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SUMMARY 


The Discovery vehicle was found to have longitudinal and Irteral aerodynamic 
characteristics similar to those of the Columbia and Challenger vehicles. The 
lateral aerodynamic characteristics of the Columbia and Challenger vehicles are reit- 
erated and the results from the Discovery flight test are added to this data base. 

The longitudinal aerodynamics resulting from the analysis of flight data from all 
three vehicles are also shown. The values of the lateral and longitudinal parameters 
are compared with the preflight data book. The lateral parameters showed the same 
trends as the data book. With the exception of Co for Mach numbers greater than 

15, ^n^^ Mach numbers greater than 2 and for Mach numbers less than 1.5, where 

the variation boundaries were not well defined, ninety percent of the extracted val- 
ues of the lateral parameters fell within the predicted variations. 

The longitudinal parameters showed more scatter, but scattered about the pre- 
flight predictions. With the exception of the Mach 1.5 to .5 region of the flight 
envelope, the preflight predictions seem a reasonable representation of the Shuttle 
aerodynamics. The models determined accounted for ninety percent of the actual 
flight time histories. 


INTRODUCTION 

During the design of the Space Shuttle, large quantities of wind-tunnel data 
wore accumulated. This information was assembled into a data base which is assumed 
to describe the aerodynamic characteristics of the Shuttle vehicle over a large range 
of flight conditions. The resulting document will be referred to as the preflight 
data book (ref. 1 ) . 

Special maneuvers were performed during the development flights of the Shuttle 
for the purpose of determining the stability and control parameters at selected 
points on the descent trajectory. Even after the Shuttle became operational, a lim- 
ited number of maneuvers were performed during the descent to the landing site. The 
responses to these maneuvers were analyzed and the stability and control parameters 
obtained were used to augment the preflight data base. The flight-derived aerodynam- 
ics were used to verify the preflight data book and to suggest areas of the preflight 
analysis that did not accurately predict the aerodynamic characteristics of the 
veh ic le • 

The results of these analyses for the lateral aerodynamic parameters for the 
Columbia and Challenger vehicles have been assembled and reports using the flight- 
determined parameters from the first eight Shuttle flights have been presented at a 
number of conferences (refs. 2-6). The main emphasis in the analysis of Shuttle 
flight data has been on the identification of the lateral aerodynamics of the vehicle 
because the preflight analyses indicated potential lateral instabilities in some 
flight regimes. The longitudinal aerodynamic parameters were determined from 
selected runs for the Columbia and Challenger Shuttle vehicles and these results will 
he presented. 

This paper will supplement the existing flight-derived parameter data base by 
adding both longitudinal and lateral results using flight test data from the Space 
Shuttle Discovery. These results will be compared with the longitudinal and lateral 
results from Columbia and Challenger. The maneuvers performed using Discovery are 
similar to those used with Columbia and Challenger. 
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acceleration measured along X, Y, and Z body a::3s, respectively, 
g units 

wing span, m (ft) 
wing mean geometric chord, m (ft) 
rolling-moment coefficient, }i^/qSh 
aerodynamic moments for trimmed flight 
pitching-moment coefficient, My/qSc 
yawing-moment coefficient, M^/qSb 
axial-force coefficient, Fx/qS 
aerodynamic force for trimmed flight 
lateral-force coefficient, Fy/qS 
trim coefficient 

normal-force coefficient, F^/qS 
vector of measurement error 

vector function representing equations of motion 
force along X, Y, and Z body axes, respectively, N(lb) 
acceleration due to gravity, 9.81 m/sec^ (32.2 ft/sec*^) 

vector function representing measurement equations 

• • 2 

moment of inertia about X, Y, and Z body axis, respectively, kg-m 
( slug-f t^) 

2 9 

product of inertia, kg-m (slug-ft ) 

cost function 

number of data points 

likelihood function 

mass, kg (slugs) 

rolling, pitching, and yawing moments, respectively, N-m (ft-lb) 
rate of roll, pitch and yaw, rad/sec or deg/sec 
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dynamic pressure, pV /2, Pa (psi) 
measurement noise covariance matrix 
wing area, m'^ f ft^) 
time, sec 

velocity along X-body axis, m/sec (ft/sec) 
input vec tor 

velocity along Y-body axis, m/sec (ft/sec) 
airspeed, m/sec (ft/sec) 

velocity along Z-body axis, m/sec (ft/sec) 

vector of states 

vector of outputs 

angle of attack, rad 

sideslip angle, rad 

aileron deflection, rad 

elevator deflection, rad or deg. 

rudder deflection, rad 

RCS control term, number of jets firing 

vector of unknown parameters 

pitch angle, rad 

roll angle, rad 

bias on roll rate, rad/sec 

quantity at i^^ time 
measured quant ity 
rotary derivatives 

static derivatives with respect to o 
static derivatives with respect to B 

control derivatives with respect to indicated quantity 
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t 

trimmed value 

X,Y,Z 

longitudinal, lateral, and vertical body axes 

Matrix exponents: 

T 

transpose of matrix 

-1 

inverse of matrix 

Mathematical 

notation: 


estimated quantity when over symbol 

• 

derivative with respect to time when over symbol 

V 

gradient operator 

Abbreviat ions 

: 

ACTP 

Aerodynamic Coefficient Tdent i f icat ion Package 

RET 

Best Estimated Trajectory 

[MU 

Inertial Measurement Unit 

LaRC 

Langley Research Center 

MMLE3 

Modified Maximum Likelihood 

PTI 

Programmed Test Input 

RCS 

Reaction Control System 

RGA,AA 

Rate Gyro Assembly, Accelerometer Assembly 

STS 

Space Transportation System 
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TEST VEHICLE 


The orbiter configuration is shown in figure 1 and key physical characteristics 
are given in Table 1. The thick, double delta wing is configured with full span ele- 
vons, comprised of two panels per side. Each elevon panel is independently actuated. 
All four panels are deflected symmetrically as an elevator for pitch control, and 
left and right elevons are deflected differentially as an aileron (6a) for roll 
control. 

The body flap is used as the primary longitudinal trim device. The elevons are 
programmed in conjunction with the body flap to follow a set schedule to provide the 
desiret! aileron effectiveness. 

The vertical tail consists of the fin and a split rudder. The rudder panels are 
deflected for yaw control and are separated to act as a speed brake to provide for 
subsonic energy modulation. The speed brake opens fully (87.2 degrees) just below 
Mach 10 and then follows a predetermined schedule until Mach 0.9 is reached. The 
rudder is not activated until Mach 5. 

Stability augmentation is provided by the aft reaction control system (RCS) 
jets, with the forward jets reserved for on-orbit attitude control and for abort 
maneuvers. The aft yaw jets are active until Mach 1, while the pitch and roll jets 
are terminated at a pressure of 20 and 10 pounds per square foot, respectively. 
Additional details of the Shuttle vehicle and its systems are given in reference 1. 


MANEUVERS 

Specially designed maneuvers were performed to obtain data for use in extracting 
aerodynamic parameters. These maneuvers were performed to obtain data at specific 
points during the descent trajectory. The test points were chosen so that 
le rodynam ic parameters could be determined along the descent trajectory to verify the 
aerodynamic model obtained from the wind tunnel tests. This verification procedure 
will add confidence to the assumed aerodynamics of the Shuttle where there is agree- 
ment and will point to areas of potential inaccuracy where there is no agreement. 

The actual forms of the inputs to be performed were developed using a Shuttle 
simulation to generate responses for various inputs and then extracting parameters 
from these responses. The control inputs that gave the best definition of the param- 
..^ters of interest were then used for the flight tests. In spite of the care taken to 
design effective inputs, because the automatic control system was active, the con- 
trols were coupled and the resulting responses were reduced in magnitude and corre- 
lated with each other and the control inputs. This Led to identif Lability and 
correlation of parameter problems during the extraction process. Additional details 
oil t)ie man'-Mjver design are given in reference 7. 


INSTRUMENTATION AND DATA PROCESSING 

As a development vehicle, the Shuttle is fully instrumented and has a number of 
redundant systems for measuring various vehicle states and controls. The instrument 
packages will be mentioned specifically. First is the Aerodynamic Coefficient Iden- 
tification Packag<=" (AC IP), an instrumentation package specifically designed to mea- 
sure rates, accelerations, and control surface positions required for parameter 
identification. The ACIP data were recorded at 172 samples per second. Second is 
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the instrumentation for the flight guidance and control system, the Rate Gyro Assem- 
bly and Accelerometer Assembly (RGA,AA), which is a source for acceleration and rate 
measurements. The RGA,AA data were recorded at 25 samples per second but are very 
noisy. The third source of flight measurements is the navigation instrumentation, 
the Inertial Measurement Unit (IMU). The IMU measurements are high fidelity but are 
recorded at only one sample per second which limits their usefulness. 

The ACIP data were the primary source for the linear and angular accelerations, 
angular rates and control surface deflections. The RCS chamber pressures were used 
to determine the jet thrust, and these measurements came from the vehicle operational 
instrumentation. 

The data considered most reliable were used to generate a best estimated trajec- 
tory (BET) for the Shuttle vehicle. The data written to tapes for the parameter 
extraction consisted of only those maneuvers considered appropriate for extraction. 
The linear and angular rates and control surface deflections came from the ACIP in- 
strumentation. The BET angular rates and linear accelerations at the start of a 
maneuver were taken as initial conditions, and the rates and accelerations were inte- 
grated over time to obtain angular positions and vehicle velocities. The velocities 
were then corrected for the effect of winds, and the resulting components were used 
to calculate the vehicle total velocity, angle of attack, and angle-of-sideslip. 

This combined data set is recorded at 25 samples per second and comprises the data 
contained on the tape to be processed by the parameter extraction software. Addi- 
tional details on the instrumentation and data processing can be found in refer- 
ences 8, 9, and 10. 


EQUATIONS OF MOTION 

The lateral-direction equations of motion used in this study are based on per- 
turbations about trimmed flight conditions and are written relative to the body axes 
shown in figure 1. The equations are 
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The results of this study are based on maneuvers performed at velocities of 
Mach I ard higher. For this reason the terms containing velocity are sufficiently 
small that tlie equations of motion are considered essentially insensitive to the ro- 
tary rlerivatives and to Cq, and ; therefore, these derivatives are fixed at 

/.ero throughout this study. 

The equations used to describe the longitudinal motions were 
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MAXIMUM LIKELIHOOD ESTIMATION 

Stability and control derivatives were extracted using the maximum likelihood 
estimator. Among other statistical properties, the maximum likelihood estimator is 
asymptotically efficient and asymptotically unbiased. This estimator consists of 
maximizing the likelihood function of the measurement errors, which is the product of 
the probability density functions evaluated at each measurement time. This approach 
requires that the form of the measurement error distribution is known; it is normally 
assumed this distribution is Gaussian. 


It is assumed the actual system can be modeled as 


X(t) - F(X,U,0,t) (20) 

Y(tJ = G(X,U,0,t.) + e., i = 1,2, ...,k (21) 


where equation (20) is a vector representation of equations (1) to (4) or equa- 
tions (8) to (13) and equation (21) is a vector representation of the measurements. 
In these equations, X is the state vector, U the vector of controls, 0 the 
vector of stability and control derivatives, t is time, and e^ the vector of 
njeasurement noise for the measurements at time 
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If it is assumed that the measurement noise is Gaussian, then the likelihood 
function (ref. 11) is 


L(Y,0) = exp 


(- J .5, lV‘i> - ""‘(V'i’ - 


( 22 ) 


where the subscript M denotes actual measurements and R is the measurement co- 
variance matrix. Taking the natural logarithm of equation (22) and multiplying by -1 
yields the cost function 


[(0) = -log L(Y,0) = 1 I - Y(t.)]'^ R"^[Y^(tJ - Y(t.)] 

i=l 


+ Y log R + 2N log 27T 


(23) 


Maximizat Lon of equation (22) with respect to 0 is equivalent to minimization of 
equation (23) with respect to 0. The last term on the right is constant relative to 
0 and can be neglected; if R is known, the second term can also be neglected for 
the same reason. Minimization of the remaining term results in solving VJ a = 0 
whichgivestheestimates 


Vl = VJ(0j), j = 0,1,2,... (24) 

Since a sequence of estimates, 0-, are obtained iteratively, the process must begin 

J A 

witli initial parameter estimates, 0^ . 

If K is unknown in equation (23), direct minimization of J( 0) with respect 
to 0 and R is complicated by the fact that R is an implicit function of 0. A 
simpler approach is to minimize with respect to 0 and R independently. Minimiza- 
tion of equation (23) with respect to R yields 


TT I - '^(t.)][Y„(t.) - Y(t.)] 

1 ■ 1 


(25) 


The procedure used here is, first, assuming R is diagonal with initial esti- 
mates for the diagonal elements, iterate equation (24) several times. Then, on each 
succeeding iteration, first estimate R using the most recent ^value of 0 in equa- 
tions (21) and (25), and then apply equation (24) once using R in J(0). This two 
step process is repeated each iteration to convergence. 
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The data from flight tests for the three Shuttle vehicles were analyzed using 
two Maximum Likelihood parameter extraction computer programs (refs. 11 and 12). The 
two programs were used independently by two investigators, and when the same data 
were processed, similar parameter values were determined. The parameter values 
resulting from processing the different Shuttle flights were plotted versus Mach 
number, and these plots form the basis for the discussion. 


ANALYSIS AND RESULTS 


In this section the results obtained in this study are discussed. These results 
are based on extracting the stability and control derivatives from 16 PTI maneuvers 
on the three flights. The time span for the measurements obtained during the maneu- 
vers ranged from 4 to 15 seconds with the measurements sampled 25 times a second. 

The estimation approach taken here is based on information contained in measured 
accelerations and rates, various trajectory parameters and the measured atmosphere. 
The method of analyzing atmospheric measurements which accounts for spatial, diurnal, 
and semidiurnal corrections is described by Price (ref. 13). This atmospheric 
in format ion is combined with onboard measurements of accelerations and rates in order 
to construct the trajectory (ref. 14) which is used for estimating the stability and 
control derivatives. 


In the results presented, lateral moment derivatives are relative to the flight 

center of gravity and were estimated with rotary derivatives fixed at zero and Cy 

. . . . ^ 

fixed at the data book value of 0.00042 per degree. The longitudinal derivatives 


'6e , 


and 


^mq fixed at data book values. All mass properties and center of 


gravity information were supplied by NASA Johnson Space Center and are shown in 

Table I. The lateral time histories fit during the estimation process were: B, 

p, r, and ay and the longitudinal time histories fit were V, a, q, and 

^or the lateral runs using the estimation method of reference 11 the weighting 

matrix (inverse of the measurement noise covariance matrix, R) was initially set to 

a diagonal matrix with the values 796.3, 234.8, 4324, 237.5, and 21820. These values 

correspond, respectively, to the measured variables B, p, r, <|), and a . Esti- 

y 

mation of R using equation (25) began on iteration 4 for each maneuver; from 8 to 
20 iterations were required for convergence. When the estimation method of refer- 
ence 12 was used, the weighting matrix was calculated based on the fit to the data 
from the initial guess at the parameter values used in the mathematical model to be 
estimated. 


The extracted stability and control derivatives will be presented in figures as 
functions of Mach number. Both flight-extracted and predicted values along with 
variations associated with the predicted values will be shown. For example, figure 2 
sliows rolling moment due to sideslip angle as a function of Mach number with the pre- 
dicted values (?) and variations (V) indicated by solid lines and the extracted 
values by various symbols. The predicted values are based on data book values, which 
are the result of numerous preflight tests of Shuttle aerodynamics (ref. 1). The 
actucal procedure for developing the predicted derivative time histories is given in 
reference 9. The variations reflect uncertainties in the data book values; they are 
based on differences between flight and predicted results for previously researched 
aircraft and extrapolated to the Shuttle configuration. The uncertainties were not 
developed for the longitudinal motion derivatives. 
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The iliscussLon of the results will be split into two parts: one for presen- 

tation of the lateral maneuver results and one for the presentation of the longitudi- 
nal maneuver results. 


at era I Maneuvers 


The results obtained from analysis of the lateral maneuvers will be examined in 
three parts: the sideslip parameters, the aileron effectiveness parameters, and the 

rudder effectiveness parameters. 

Side-Slip Parameters . Cyg - The side force due to sideslip parameter is shown 
plotted versus Mach number in figure 2. This parameter tended to be difficult to 
identify, and the Cramer-Rao bounds were generally greater than 10 percent of the 
extracted value for the parameter. If there is a trend, the Challenger tends to 
show less At the lower Mach number o becomes more negative, but how much 

is not clear. 


The trends for the Discovery vehicle seem to follow the preflight predictions 
well except for the scatter at the lowest Mach numbers. The results show no signif- 
icant differences between the Discovery vehicle and the Columbia and Challenger vehi- 


cles (fig. 2). For the Cyg 
flescription of this derivative 


parameter, the preflight prediction offers a reasonable 


Cj^g - The rolling moment due to sideslip parameter was well determined, and the 

Cramei — Rao bounils were generally less than 2 percent of the extracted parameter 
value. The values of Co- were less negative than predicted at the highest Mach 

P ^ 

numbers (fig. 2). The values determined from the Discovery flight data consistently 
followed this trend established by the other vehicles. In fact, the Discovery re- 
sults were generally more consistent run-to-run than those of the other vehicles for 
most of the Mach range. Below Mach 5, there is considerable scatter in the C^ 

values extracted for all the vehicles (fig. 2). In this Mach range, the generaf 

Cft value near zero for 


trend agrees with the preflight predictions and indicates a ^ 


e 


Macli number less than 1. The preflight predictions are a reasonable representation 
of Cq for the mid to low Mach numbers (fig. 2). 

- The yawing moment witii sideslip angle parameter is shown plotted versus 
M/ic.li muilher in figure 2. The flight determined C^ generally agrees with the pre- 
i light predictions for all vehicles. Only in the Mach 5 to Mach 10 range do the 
values tend to be m()re negative than predicted. The trends below Mach 5 seem to be 
more clearly defined for all vehicles than for most of the other parameters extracted 
(li;> 2). Again, the C_ values extracted from the Discovery flight data agree 
* ^ 3 

very well with the values determined from the flight data of the other two Shuttle 
veh i c les . 


Aileron Control Effectiveness. C 




The rolling moment due to aileron versus 


Macli number is shown in figure 3a. The values extracted from flight data for all 
three vehicles showed the same trends as the preflight predictions. While the 


Cramer-Kac bounds for C, 


"6a 


were, in general, less than 2 percent of the extracted 


values, the scatter in the extracted values implies that the parameter may not be as 


well determined as 




and (fig* 2). The Discovery vehicle aileron showed 
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less effectiveness than the other vehicles at high Mach numbers, but the limited data 
make conclusions difficult. As a general conclusion, for all Mach numbers, the 

preflight estimates seem to be a good representation of C« . 

^6a 

- The yawing moment due to aileron parameter plotted versus Mach number is 

* . . . 

sliown in figure 3b. The parameter was not especially well determined at the high 

Mach numbers, and the Craraer-Rao bounds are on the order of 10 percent of the 
extracted parameter value. In the Mach 18 to Mach 5 range, the determination was 
better, but the Cramer-Rao bounds were still high. Below Mach 5 showed the 

same scatter that has been seen from most of the parameters. The values of ^ng^ 

for tlie Discovery vehicle generally followed the same trends as those seen for the 
Columbia and Challenger vehicles; however, they seem to indicate slightly less aile- 
ron effect at the higher Mach numbers (fig. 3b). The values determined from flight 
data for all three vehicles followed the same trends as the preflight predictions. 
Hence, these predictions are a reasonable representation of the value of 


Rudder Control Effectiveness . rolling moment due to rudder parameter 

plotted versus Mach numbers is shown in figure 4a. The parameter was well determined 
and the values extracted from flight data for each of the Shuttle vehicles showed 
trends similar to those of the preflight predictions. The values determined 

using the Discovery flight data agreed well with the preflight predictions, even at 
the lower Mach numbers. Again, the parameters determined using flight test data are 
similar in value and trend to the preflight predictions indicating that the latter 

are reasonable for the parameter. 

^6r 

C^g^ - The yawing moment due to rudder parameter plotted versus Mach number is 

shown in figure 4b. The parameter was well determined and did not show much scatter. 
The values determined using the Columbia and Challenger vehicles showed similar 
trends and values for ^ngj.* values determined using the Discovery flight data 

agree well with those obtained using flight data from the other vehicles. However, 
for all vehicles C^ was consistently less effective than predicted above Mach 2. 


In summary, the parameter values determined using data from the three vehicles 
agree very well. In most cases, the parameters, with Cramer-Rao bounds approximately 
?. percent of the extracted value, agree well with the preflight predictions. As a 
further check, several places were found on the Discovery flight tape where there 
were aileron trim settings prior to maneuvers. For each aileron trim there was a 
trim sideslip angle and the ratio gave an indication of the C^ to Co ratio 

:3 . ^6a . 

that might be expected from the maneuver. Since there were trim points at different 
Mach numbers, any variation in this ratio with Mach numbers should be indicated. The 

trim ratios, estimated Co to Cn ratios and the predicted C« to C- 

^6a . ^6a 

ratios, are all given in Table II for several Mach numbers where the trim was well 

enough established to make the calculations. The data in the table show that the 
parameter ratios calculated from trim were in fair agreement with ratios determined 


using the extracted parameters, and both showed reasonable agreement with the pre- 
flight prediction ratios. The trends with Mach number were also predicted by the 
trim settings, indicating that the trends determined from the flight test for 
and ^£53 seem reasonable. 
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Longitudinal Maneuvers 


Attempts to determine longitudinal aerodynamic parameters from flight data for 
the Shuttle vehicle have been only partially successful. Several factors contribute 
to this limited success. First, the perturbations that are allowed by safety consi- 
derations are so small that the vehicle cannot be disturbed sufficiently to have 
responses that are large enough for good parameter Identification. Secondly, only 
the control surfaces are available to perturb the vehicle and these move so slowly 
that true step inputs cannot be used to disturb the vehicle. Responses to the slower 
Inputs can lead to correlation problems during parameter extraction. Hence, a 
limited discussion of selected longitudinal parameters will be presented. 

Because of Ident I f iabi 1 ity and correlation problems, the longitudinal parameters 
will first be discussed as a group and then individual parameter characteristics will 
be discussed. The discussion will center on three longitudinal parameters; 

the change In force along the **Z’* body axis with angle-of-attack change; the 

pitching moment due to angle-of-attack change; and ^m^e^ pitching moment due to 

eleven change. The motion of the Shuttle about some trim conditions resulting from a 
control Input can he described using a set of equations containing (trim 

coefficient) and the three parameters under discussion. This reduced mathematical 
model will describe about 90 percent of the motion about the trim condition. Since 
the parameters and describe the majority of the vehicle motion, 

they can be Identified with the greatest confidence, and only these parameters will 
he discussed In the paper. The three parameters are plotted versus Mach number and 
are shown as figure 5. 

The ^m^^ parameter was best defined with a Cramer-Rao bound less than 3 per- 
cent of the extracted parameter value. Cramer-Rao bounds that were about 

^ percent of the extracted value while had bounds that were about 10 percent of 

the extracted value. C and Cm were correlated In almost all cases. In spite 

of the correlation and poor Ident i f lab i 1 I ty problems previously mentioned, 

, and C generally showed trends similar to those of the preflight data book 
?e "’a 

(fig. 5b and 5c ) . 

“■ The parameter values determined showed that above Mach 15, was 

);eneralLy more negative than predicted (solid curve) when the Maximum Likelihood 
estimator was used (fig. 5a). In many cases below Mach 15 the values obtained using 
t}ie Columbia flight data tended to be less negative than predicted while the values 
obtained using the Challenger data tended to be more negative than predicted. The 
trends seen using Columbia data were also evident in the results of reference 15. 
Values of calculated from steady-state relations showed trends close to those 

of tlie pref light predictions (fig. 5a). Most of these calculations were made using 
data from Columbia flights. 

The values for C^^ obtained using the data from the first Discovery flight are 

similar to those for the Challenger vehicle. However, when the flight data from the 
second Discovery flight were used to determine C^^ above Mach 5, the values ob- 
tained were consistently less negative than the values obtained using the data from 
the first Discovery flight and from the Challenger flights (fig. 5a). An examination 
of the different sets of data showed that pulse-type commands were used to excite the 
vetiicle for Longitudinal parameter extraction during the Challenger flights and also 
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for the first Discovery flight* During the second Discovery flight doublet commands 
were used above Mach 5 to excite the vehicle. The actual eleven surface motion was 
not a true doublet because the automatic control system modified the input. A sample 
time history of the elevon inputs for the first and second Discovery flights is shown 
In figure 6. The parameter values extracted for each input are given in Table II. 

The value is clearly significantly different. In most cases the Cramer-Rao 

bounds were about 5 percent of the extracted parameter value implying some confidence 
In the values extracted. 


When was fixed at the preflight predictions, ^ 5 ^ were unchanged or were 

close to the pref light values In a majority of the cases examined (fig* 5a). This 
new set of parameters still described about 90 percent of the perturbation from trim 
motion of the Shuttle. This implies that the restricted inputs could not excite the 
vehicle enough to Identify all of the parameters in the assumed mathematical model, 
and that although seemed reasonably well identified, a change of over 50 per- 

cent in the parameter value did not have a serious effect on the fit to the data or 
the values of ^mg^ fixing for selected runs throughout the Mach 

range, a model closer to the preflight predictions that still fits the flight data 
could be obtained. Below Mach 5, however, some questions still remain as to the best 
model to represent the Shuttle longitudinal aerodynamics and this will be addressed 
In the section discussing C and C 


and 


m 


6e 


- Since these parameters are correlated for most of the Shuttle 


runs examined, they will be discussed together. At the higher Mach numbers there was 
considerable scatter In the extracted values for C , but the trends generally 

followed those of the preflight predictions (fig. 5b). Below Mach 5, the values 
extracted for C tended to show less elevon effectiveness than predicted. For 

the Columbia vehicle this trend was also shown in reference 15. However, in the Mach 
greater than 5 range, when was fixed at the preflight predictions, 

. ^ 

generally changed very little. While C variations were less that 10% for most 


a 


cases there were some significant variations seen. This was not unexpected because 
C was not well Identified. Since is considered better defined than the 


m 


6e 


otlier two parameters, the fit to the data did not suffer significantly. 


For Mach less than 5, the trends of the extracted values of 


6e 


and C„, 


followed the trends of the prefllght predictions down to Mach numbers less than 2 
(fig. 5). Below Mach 2 the prefllght predictions show large variations in and 

(fig* ^)* The magnitude of the values generally followed the trends of 

the prefllght estimates. However, the extracted values of C 


""6e 


were generally 30 


to 50 percent less than the predicted values. In this Mach range fixing did 

not significantly affect In most cases, but some significant changes were seen 

in C (fig* 5b). The most serious difference between the preflight predictions of 

^ ot 

the longitudinal parameters and the flight determined parameters is the reduced 

elevon effectiveness as indicated by C around Mach 1. This is illustrated by 

^ 6e 

considering the situation where both and are set to the prefllght pre- 

dictions, then the values of C extracted were as much as 100 percent different 

from the preflight predictions, and in many cases, became positive. The fit also de- 
graded for these cases. Although there are only five runs near Mach 1 where the 
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predLcted greatest magnitude, the preflight predictions imply more 

control effectiveness than was found using the flight test results. 


Above Mach 5 the Discovery vehicle results show the same trends for 


and 


m 

xc i te the veluc Le • 


as the Challenger vehicle when the same control inputs were used to 


[» 


C 7 was consistently more negative than the preflight pre- 


<lictions,. However, when a different input form was used, the values were consi- 
derably less negative than predicted. For the Columbia and Discovery maneuvers where 

Cy was less negative than predicted, C and C still generally showed the 
'"a , . 6 e 

same trends as seen with the Challenger and first Discovery flights. 


In all cases using Discovery flight data taken above Mach 5, fixing 

preflight values resulted in a mathematical model that described the motion of the 
vehicle and that showed trends similar to those of the preflight predictions 
(fig. 5). 


For C,„. and C^ in the Mach number less than 5 range, the values estimated 
‘■‘oe "‘CX 

using Discovery data were close to the preflight predictions except around Mach 1, 
and fixing did not change these values significantly. As previously discussed, 

the values f)f C,„ . extracted from flight data in the Mach 1.5 to 0.6 range, were 
oe 

significantly different from the preflight predictions. However, for the rest of the 
Mac)» range, the preflight predictions should be considered a reasonable representa- 
tion of the aerodynamics, not only of the Discovery vehicle, but for the Columbia and 
Challenger vehicles also. 


CONCLUSIONS 

The Discovery, Challenger, and Columbia all exhibited similar lateral and longi- 
tudinal aerodynamic characteristics. The preflight predictions are a reasonable 
representation of the Shuttle aerodynamics for all vehicles except for three cases. 
The first is for in the Mach number greater than 15 range where the flight- 

determined parameter value is smaller than predicted. The second exception is that 
for Mach numbers greater than 2 was determined to indicate less rudder effec- 
tiveness tlian predicted. Thirdly, the. Mach 1.5 to 0.6 range indicated 

.il)proxiina te ly thirty percent less eleven effectiveness than predicted. 
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Table 1 Comparison of Trim Curve and 
Parameter Ratios 


Mach 

Number 


Aileron Trim/ 
Beta Trim 


Flight Predicted 


'Z 


3 


^3 


'6a 


'6a 


1 

.93 

.91 

.93 

3 

2.32 

2.04 

1.96 

4 

2.5 

3.23 

2.0 

6 

1.82 

1.67 

1.43 

22 

1.41 

1.39 

1.32 

Table II 

Parameters 

Extracted from 

Di f ferent 


Discovery Inputs at Mach = 8 
for Inputs of Figure 5 



Parameter 

Input 1 

Input 2 



-1 .03 

-.96 

0 




00 

. 

. 

00 

a 




-.11 

-.15 

a 




-.16 

-.13 

6e 
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Senses indicated 
are positive. 


Z 

w 

w 


Figure 1 



- Schematic of Shuttle vehicle showing body axes and 
positive senses of accelerations, rates, velocities, 
moments and angles. 







O Discovery 
O Coluabl* 

0 Challenger 

A Steady State Calculetlona 
Shaded eyabola Indicate 

1 I ^ predicted valves a 



(a) Change in nonal force with alpha 






(c) Elevon effectiveness 

Figure 5.- Dominant longitudinal derivatives affecting the abort period 
■otlon of the Shuttle versus Mach nunber^( derivatives per degree). 
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